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ABSTRACT 

Variations  on  the  traditional  cost  minimization  of  continuous  review 
formulation  are  investigated  in  an  effort  to  improve  service  as  measured 
in  terms  of  time-weighted  shortages  per  unit  time.   It  is  proposed  that 
the  minimization  of  time-weighted  shortages  per  unit  time  will  improve 
service  in  current  Navy  Supply  Operations.   Various  models  are  presented, 
without  reliance  upon  unknown  parameters  such  as  order  cost  and  carrying 
cost,  with  necessary  conditions  and  solution  algorithms. 
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I.   INTRODUCTION 

Inventories  exist  to  provide  service  to  customers  by  satisfying  their 
demands  from  on-hand  material.   It  follows,  then,  that  a  reasonable 
objective  of  inventory  management  is  the  maximization  of  service  provided 
which  is  achieved  by  minimizing  stockouts.   In  particular,  total  time- 
weighted  shortages  is  thought  to  be  the  desired  objective. 

In  pursuing  this  objective,  the  manager  of  a  realistically  large, 
multi-item  inventory  system  has  a  number  of  constraints  imposed  on  his 
"when  to  buy  and  how  much  to  buy"  decisions.   The  stock  points  of  the 
Navy  Supply  System  have  investment  and  reorder  workload  constraints  which 
are  real  and  binding. 

The  classic  variable  cost  minimization  formulation  is  the  most  used 
method  for  solving  this  inventory  problem.   Multi-item  problems  are 
usually  solved  by  assuming  that  they  can  be  dealt  with  as  a  series  of 
independent  single  item  problem.   In  the  presence  of  binding  constraints 
on  a  population  of  items  this  approach  is  not  applicable.  Additionally 
the  cost  minimization  formulation  requires  the  estimation  of  cost  para- 
meters which  are  arbitrary  or  at  least  very  difficult  to  estimate. 

As  a  consequence  of  this  argument,  a  series  of  models  are  formulated 
for  multi-item  policies  subject  to  investment  and  reorder  workload  constraints. 
These  models  do  not  employ  the  standard  ordering,  shortage  costs.   This 
approach  was  suggested  by  A.  P.  Tully  [_lj  . 

In  the  next  section,  the  problem  formulation  and  the  general  models 
are  developed.   Section  III  develops  the  single  item  model  as  preparatory 
to  studying  multi-item  case.   Section  IV  presents  a  simplified  multi-item 


formulation  in  whicn  only  the  items  reorder  points  are  decision  variables. 
The  general  multi-item  continuous  review  model  is  developed  in  Section  V. 
In  the  final  section,  we  summarize  and  state  some  tentative  conclusions. 
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II.   FORMULATION 

It  is  desired  to  formulate  inventory  decision  rules  for  multi-item 
inventories  subject  to  specific  constraints.   The  inventory  decision  rules 
will  be  of  the  reorder  point  -  reorder  quantity,  continuous  review  type. 

As  suggested  by  the  introduction,  the  formulation  to  be  used  involves 
the  minimization  of  total  time-weighted  shortages  subject  to: 

(1)  Total  average  investments  costs  /L    investment  limit;  and 

(2)  total  number  of  orders  <d  reorder  workload  constraints. 

Note  that  such  a  formulation  would  not  be  based  on  minimizing  variable 
costs. 

The  specific  form  of  the  model  depends  upon  the  assumption  about  the  item 
demand  characteristics  and  expressions  for  the  total  average  on-hand  inventory 
level  and  total  number  of  buys  per  unit  time.   The  first  assumption  is  the 
distribution  of  lead  time  demand  is  normal  (ife,  »i  )  for  all  items. 

The  following  notation  is  used  throughout  the  paper.   For  the  i-th  item  let; 

Ci  =  item  unit  cost  in  dollars; 
fal   =  mean  demand  per  unit  time  in  units; 
M.i   =  mean  lead  time  demand  in  units; 

Oi  =  Standard  deviation  of  lead  time  demand  in  units; 

$(r^)  =  probability  that  lead  time  demand  exceeds  r; 

r^  s  reorder  point;  and 

Qi  =  reorder  quantity. 
Also  let; 

K,  =  investment  limit;  and 

Kj  =  reorder  workload  constraint. 
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With  a  continuous  review  inventory  policy  an  order  is  placed  after 
the  demand  of  Q  units.   It  follows  then  that  the  average  number  of  orders 
per  unit  time  is  -^—   .   For  a  multi-item  inventory  with  N  items,  the  total 
expected  number  of  orders  placed  per  unit  time  is 


Z%i  (II. 1) 

Qi  . 
i=l 

Total  inventory  investment  is  the  priced-out  value  of  the  total 

expected  on-hand  inventory.   As  shown  by  Hadley  and  Whitin  \_2A     ^or 

continuous  review  the  expected  on- hand  quantity,  E(OH)  is  given  by 

E(OH)  =  r  +5  .fj,   +   B  (Q,r), 

where  B(Q,  r)  is  the  expression  for  the  expected  shortages  at  any  point 
of  time.   If  lead  time  demand  is  normally  distributed  it  can  be  shown 

OQ   that 

B(Q,  r)  =~  [f(r)  -  B  (r  +  Q)  ]  ,  (II. 2) 

where       B(r)  =  I  &   +  (r  -A  )2  1  $  (  S^>  f  (^->^> 

1„2 

e 

f2Tf 

$(r)  =  \      <j?(X)dX. 


(()(X)  =  ~L  e"  2X 


The  expected  on-hand  quantity  expression  can  be  simplified  by  omitting 
the  B(Q,r)  term,  and  this  approximation  is  reasonable  if  the  risk  of  stock 
out  is  not  too  large.   This  assumption  is  employed  throughout  the  thesis. 
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With  this  assumption,  the  total  inventory  investment  is  then  given  by  the 
expression 


N  n. 

XT"  Ci  (  ri  +  ^  -  ytti).  (II. 3) 

i=l 


The  expected  number  of  backorders  at  any  time  may  be  explicitly 
determined  from  the  steady  state  probability  distribution  for  negative 
net  inventory  levels.  Hadley  and  Whitin  [_2_J  used  this  approach  and 
showed  that  when  lead  time  demand  is  normal,  the  time-weighted  shortages 
expression  is  given  by  Eq  (II. 2). 

If  the  risk  of  stockout  is  small,  then  the  expression  for  the  time- 
weighted  shortages  can  be  simplified  by  ignoring  the  B(r+Q)  term,  which 
yields  expected  time-weighted  units  short  per  unit  time  for  the  ith  item 
as: 

-^-  p(ri)  (II. 4) 

The  objective  can  now  be  stated  as  the  minimization  of  the  time- 
weighted  shortages  for  the  entire  inventory,  and  the  formulation  is: 


.N 
min  Z  = 


I    iff 

i=l 


(ri)  (II. 5) 


subject  to: 

N 

I 

i=l 

Ci(ri+^  -  jUi)   £    Kx  , 

N 

I 

i=l 

Qi       -   2  » 

Qi> 

0,  and  r^   unrestricted. 
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(II. 6) 


(II. 7) 


Note  that  the  investment  constraint  will  always  be  active,  given  the 
objective  function  used,  but  that  reorder  constraint  may  or  may  not  be 
active  in  a  given  problem. 
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III.   SINGLE  ITEM  B(Q,r)  MODEL 


The  basic  continuous  review  formulation  for  one  item  with  normal 

2 

(U,  d   )  lead  time  demand  can  be  written  as 


min: 

-T    p(r) 

subject   to: 

C   (r+  |  -JUL) 

4iH  . 

Q 

z  k2  , 

(III.l) 

(111. 2) 

(111. 3) 


Q  >  0,  and  r  unrestricted. 

A.   NECESSARY  CONDITIONS 

To  solve  the  single  item  continuous  review  model,  set  up  the 
Lagrangean  function 

L(Q,r,4t,a)     =    \     £(*)   +    fl[c<r+S     -  A)-*±] +8l3[' fyl     (m.4) 

Taking  the  partial  derivatives    ~r  ,  —  ,  r-r  ,  and  =tk   , 
and  setting  the  resulting  expressions  equal  to  zero  yield; 

3jj--±jP  (r)+^-i|  =  0    ,  (III.5) 

|?"  ^CCt-)i^>)-*i^)D  +41C-0,  (III.6) 


|£  =  C(  r  +  |  -jtt  )  -  K^  -  0  ,  and  (III. 7) 

If  =t-k2  =  o  •  <m-8) 
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These  equations  may  be  rewritten  as 

Q2  ,  2  Ep(r)  +  6X3  ,  ttI>9) 

<t<r^r)   -  (r-A)$  (^)  =UCQ,  (in. 10) 

C  (r+|  -JUL)   =  Klf  (III. 11) 

and    ^L  =  K2  .  (III. 12) 

These  are  the  necessary  conditions  for  solution. 

B.   ITERATIVE  SCHEME 

If  the  reorder  constraint  is  active,  the  reorder  quantity  is  determined 
from  the  equation  (HI. 12) 

£  =  K     which  implies: 
Q    2 

Q=4i  (III. 13) 

2  * 

Note  that  in  this  case,  equation  (III. 11)  can  be  solved  for  r  yielding 
r  =-^  -~+  A   •  (III. 14) 

Hence  (Q*,  r*)  are  uniquely  determined  from  equations  (III. 13)  and  (ill. 14). 
If  the  reorder  constraint  is  not  active,  the  reorder  quantity  is 

determined  along  the  line  which  is  equation  (III. 11).   It  is  observed  that 

Kl   if 

the  boundary  conditions  are  as  follows:   when  r  =  0,  Q  =  2  ( —  +/**-),  and 

M 
K1 

when  Q  =  0,  r  =  — -  +/1.      Plotting  these  conditions,  Figure  1  is  obtained. 
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r  = 


(rmax,  Q  =  0) 


(r  max  if  reorder  constraint 
active) 


(r+§./t)  =  K 


Figure  1.   Boundary  Conditions  for  Single  Item  Case. 
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it  if  -3f_  i 


Kl 


It  is  realized  that  if  -^-  is  greater  than  2(—  +A  )>  then  the 


K 


2 


C 


investment  and  reorder  constraints  are  mutually  inconsistent  (infeasible) . 

0  1 

If  all  solutions  lie  on  the  line  C(r  +  ■*  -JUL   )  =  K  for  Q  >  ~  and 


K 


^  3  +/i- 


4. 


C     2K   ,  then  pick  up  a  Q,  solve  for   from  the  equation  (III. 11), 
evaluate  v'ry    and  continue  to  iterate  on  Q  until   ^ -*   minimum  is 

Q  Q 

found.   Here  parameterizing  K  will  give  shortages  as  a  function  of  the 

investment  limit. 

An  alternative  approach  would  be  to  use  a  Lagrangean  function  with  the 

boundary  conditions  \  £■  JL   Q  L  2( _i  +/1  )  I   and  the  equation  (III. 9) 

LK2  "   '     C   ^   J 

with  9=0.   From  these  equations,  it  is  determined  that 


%  9  Ki      2 

2(-~+/t)  C  (III. 15) 

Hence  a  double  iteration  is  required  as  the  function  of  f^  ,  Q,  and  r  ;  i.e., 
iterate  Q, r  for  a  given  value  of  the  multiplier.   The  solution  is  the 
minimum  time-weighted  shortages  for  the  multiplier  used;  minimum  for  the 
investment  level  (K. )  imputed  to  that  multiplier  value  4y  .   Then  a  sug- 
gested procedure  would  be  a  binary  search. 

Plotting  shortages  vs.  multipliers  on  the  Q, r  and  shortages  planes, 
the  different  level  of  shortages  will  be  obtained. 

C.   EXAMPLE  OF  THE  SINGLE  ITEM  MODEL 

Consider  an  item  with  its  distribution  of  lead  time  demand  normal 
(/*1  <?2).   Let  fl*  =   125,  d2   =  25,  %  =  500,  C  =  12,  K  =  720,  and  K2  «  8. 
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The  problem  is 


minimize 


Z=  ±  p(r) 
bject  to:   C(r+  |  -/t)  X  K  , 


*  £        K 

o    *     V 


QVO,  and   r  unrestricted. 
From    -  /.   K  ,  it  is  determined  that  Q  %  —^  =  62.5.   Hence  if  the 


8 


reorder  constraint  is  active,  Q*  =62.5,  r  is  obtained  from  (III. 11)  as 
r*  =  144,  and  Z*  =  .606. 

If  the  reorder  constraint  is  not  active,  then  a  double  iteration  is 
required  as  the  function/^  ,  Q,  and  r  .   The  solution  is  the  minimum  time- 
weighted  shortages  for  the  multiplier  used.   Iterating  on  Q,  and  ignoring 
the  reorder  constraint,  produces  the  following  results: 


Time-weighted 

Expected  number 

shortages  per 

of  units  short 

Q 

r 

unit  time 

per  unit  time 

10 

170 

0.3198 

17.8446 

18 

166 

0.027 

14.5262 

36 

157 

0.333 

16.2599 

50 

150 

0.471 

20.8288 

62 

144 

0.606 

25.0432 

76 

137 

0.904 

32.9660 

90 

130 

1.247 

42.6242 

105 

122 

1.741 

53.6828 

Table  1.   Time-weighted  Shortages  and  Expected  Number  of  Units 
Short  for  Single  Item  Case. 

Note:   Two  approaches  seem  to  follow  together. 
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IV.   SIMPLIFIED  MULTI-ITEM  B(Q,r  )  MODEL 

The  basic  continuous  review  formulation  for  the  multi-item  problem 
was  given  as: 


minimize    Z  =  ST~     ft(ri)  (IV.  1) 

L.  Qi 

i=l 


subject  to:    N 


V"  Ci(  ri  +  ^-  -  Ai)  £  K   , 


(IV. 2) 


1=1 

N 


|   ^.  <IV'3> 


1=1 

Qi  ^  0  ,  and  r^   unrestricted. 

Suppose  the  order  quantities  are  fixed  by  some  other  criterion. 
Specifically  the  assumption  is  made  that  order  quantities  are  determined 
from  the  equation   Qi  =  G  l^r  ,  (IV. 4) 

where  G  is  a  constant,  which  is  assumed  to  be  the  same  for  all  items. 

From  the  reorder  constraint  (IV. 3)  and  the  assumption  (IV. 4)  it  follows 
that 


N 


1=1 


£—  -  K2  ,  (IV.  5) 


or 


N 


J%iCi 


K2    . 

The  determination  of  G  then  fixes  the  order  quantities  from  equation  (IV. 4) 
and  eliminates  one  set  of  decision  variables  from  the  problem,  i.e., 
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N 

2—   JXiCi 

i=l 


K, 


(IV. 6) 


Substituting  (IV. 4)  into  equation  (IV. 2),  it  follows  that 


N 


i=l 


Ci  (  ri  +  -7T 


2  4  Ci 


Xi 


J   -Jttl)   Z,  K 


Equation  (IV. 7)  can  be  reduced  to  the  form 


(IV. 7) 


N  N         N 

V"   C±'*t&\-     —  ^>jXiCi  +^Ci^i  =  K' 


(IV. 8) 


1=1  1=1       1=1 

Now  the  multi-item  B(Q,  r)  Model  with  fixed  Q's  can  be  written  as 


N 


minimize  Z  = 


subject  to: 


z 

1=1 


p 


(£i) 


Qi 


N 

H     ciri  ^k;  . 
i=i  i 


and  ri        unrestricted, 

N 
where  Qi  =       /       JXiCi 

l-T  * 


2K, 


%i 


M  Ci 


=  G   \Tci 


A.   NECESSARY  CONDITIONS 

To  solve  the  simplified  model,  set  up  Lagrangean  function 


(ri'V'i  -rJS   P<rt>+4[^ci-i-rJ 


Ci  ri   -   K  (IV. 9) 
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Taking  the  partial  derivatives  with  respect  to  decision  variables  and 
setting  these  expressions  equal  to  zero  yields 

Ul-^JI  [<*■****&*>    -<M(-£f^>]   *H«  "  0,  and 

|t-  -  £  ci  rl  -  K;  -  o. 

1=1 

Thus  the  conditions  for  solution  to  the  problem  are 


or 


*  -    -^        [rfi  4t^j&)  -   (  ri-A>  $  (-^)]        ,     (Iv. 
gJciAi        u 


10) 


N 
and    V"   Ci  ri  =  k'  (IV.  11) 

i=l 


B.   NUMERICAL  SOLUTION  SCHEME 

In  general,  these  equations  cannot  be  solved  in  closed  form.  A  numerical 
solution  procedure  is  suggested.   Let  us  consider  a  numerical  method  of 
solving  equations  (IV. 10)  and  (IV. 11).   Observing  equation  (IV. 10),  which 
is 

note  that  the  right-hand  side  has  a  lower  bound  of  0  since  4|*fe0  from 
G  >0,  Ciii>0  and  ^i <j>  (-* 'pf— )  -  (ri-y^i)  $  (  ^   )   0  for  all  ri. 

Specifically fy  =  0  when  ri  becomes  infinite  and  0^i<ft(  TiJ'i)    -    (ri"Ai)» 
(p  ( •J'     )  =  0,  br   r^  =(JD  violates  the  investment  constraint. 


22 


Hence  for  the  initial  value  of  jh     it  is  reasonable  to  start  with 
r.  =  0,  then  this  implies 


or 


'i^i  L  °i  aiJ 


JLmin  -J      fl'.Af-^  +  jtt.  J  (-  —  ~)\      for  all  i  ; 


ffl   »  g  implies  that  there  is  at  least  one  r  at  its  lower  bound  of  zero. 

Hence4).=  —  is  a  convenient  starting  point.   Then  a  suggested  solution 

2  C 

procedure  would  be  to  begin  at  4L  =  — ,  solve  equation  (IV. 10)  for  the 

r  's  and  compute  the  value  of  constraint  using  equation  (IV. 12),  which  is 

N 

Y~   C^  =  H.  (IV. 12) 

i=l 

g 

A  bisection  search  will  be  used  again.   If  H^  K  ,  increase^  by  . 

If  H  ^  Kj*  ,  decrease  ^  by  -5—  .  Recompute  the  r^'s  and  the  value  of 
constraint  using  equation  (IV. 12).   If  the  increase  (or  decrease)  of  fy 

has  not  caused  the  change  of  inequality  sign,  increase  (or  decrease) /ft 

C 

by  the  same  amount  _£_  .   If  the  sign  of  inequality  has  changed,  then 

4   O 

reduce  the  increment  to  — —  and  increase  (or  decrease)  ^f  ,  solving  for 

the  r^'s  at  each  value  of  %  and  computing  the  value  of  H  until  the 

sign  of  the  inequality  switches  again.   Continue  until  H  =  K  '  or  until 

H  is  within  some  tolerable  limit  of  K,'  .   This  approach  will  converge  to 

the  optimal  solution  rapidly. 
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From  the  Kuhn- Tucker  theorem  [_3j  ,  If  we  have  a  convex  objective 
function  and  a  convex  constraint  region,  the  necessary  conditions  are  also 
sufficient.   Since  the  constraint  under  consideration  is  linear  in  r  ,  the 
region  is  convex.   To  show  Z(ri)  is  convex,  consider  the  equation  of  the 
expected  time-weighted  shortages, 


Ci    P(ri). 


Zi  =_1_ 

G  ^i 

2  a 
Now  if    — ^  0    for  all   ri,  then  Zi  is  convex.   Taking  partial 

derivatives, 


*&Zi     _ 


HI   [(  ^-AD$(-^i)  -   *i*(^i)]  <o, 


fcri  G 

and 

^ii-=      J-fcT    a,    rl-/*l.  rTvrn 

©n2  g  J;j£    2<— oq~ >     •  <IV-13> 

Equation  (IV. 13)  will  always  be  greater  than  or  equal  to  zero.   Under 
these  conditions   Zi  is  convex.   It  follows  that  Z  is  convex  since  it  is 
the  sum  of  convex  functions. 

C.   EXAMPLE  OF  THE  SIMPLIFIED  MULTI-ITEM  B(Q,r  )  MODEL 

Let  us  consider  an  inventory  of  three  items.   It  is  assumed  that  the 
distribution  of  lead  time  demand  is  normal  with  mean  M*i.   and  variance  ^i 
for  the  itn  item.   Let  the  item  data  be  as  follows: 

Item  1         Item  2         Item  3 


Ok 

1000 

1500 

2000 

Ci 

1 

10 

20 

*K 

100 

200 

300 

•i2 

100 

100 

200  ,  and  let 

K 


=   $8,000  and  K2  =  15, 
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Now  the  problem  is 

3 
minimize  Z  =  -L-   ^~  fit    3  (  ri) 


G    fe  tfV 


subject  to: 


3  3  3 

i=l  i=l        T=T 


where 

G   = 


K2 


=  23.6065. 


order  quantities  are  determined  from  equation  (IV. 4), 

Ql  =  746.5022, 

Q2  =  289.1190,  and 

Q3  =  236.0648  . 

With  these  orde-:  quantities,  the  problem  now  can  be  stated  as 


minimize  Z 

3 
i=l 

1 

rcL 

4^i 

23.6065 

subject  to: 

3 

jT     Ciri  ^  8000  +  y-  ci/ti  -  ^  ^^   =  11,920. 
i=l  At  i=l 


5078 


i=l  i=l 


Notice  that  the  solution  to  the  above  problem  is  that  vector  r  such  that 

%  -  — 1—  [<M(-J:Vfi>  -  <  ri-AO^C-^f^J   '  f°r  aU  l«  (IV'10) 
G  ^Ciri  L 
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For  each  value  of fy   ,  there  will  be  some  £  r  }  for  which  equation 
(IV. 10)  is  satisfied.   However,  having  the  convex  objective  function, 

there  exists  only  one  jrij  for  all  i  such  that  equations  (IV. 10)  and 

g 

(IV. 11)  are  satisfied  at  the  same  time.   For  the  initial  value  of  /fl   =  —  , 


us 


ing  g  =  -1-  -in  -1      U+l-   ^L>  +  *i  $  (-  £f  )1 

for  all  i  when  ri  =  0,   %o   =  0.0318  will  be  used.   Decreasing  ^  from 

0.0318, ^p  =  0.0055  is  obtained,  which  implies 

rl  =  234.3750, 

r2  =  264.0625,  and 

r3  =  453.1250. 

Checking  the  constraint,  it  is  found  that 

3 

/        Ci  ri  =  11937.5000, 

i=l 
which  is  within  0.2  percent  level  of  K,'  . 

The  expression  for  the  time-weighted  shortage  per  item  per  unit  time 

is 

The  time-weighted  shortage  per  unit  time  is  computed  as 

*  "  *Ab  JX  [  *  H  +  ^M'UO-  ¥  C240-100«)] 

\   [l002  +  (285-200)2]$(I|§)-^(285-200)(|)(T|§)J 


"L        23.6065 

JlOOO 

=  0.2179 

> 

n       =              1 

rrr 

~z       23.6065 

J  1500 

=  2.7634   , 

and 
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Zq   ■  ^ QlT!    f  2002+ (285-300)  2I  $(--H)    -^00    (285-300)<f>  (-JJ> J 

3        23.6065    J2000  L   2    I  J*      200  2  •        200;J 

=   10.5232 

Total  time-weighted  shortages  per  unit  time  for  the  entire  inventory 


are 

Z  = 

i=l 


V"  Zi  =  13.5045. 
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V.   GENERAL  MULTI-ITEM  B  (Q,r  )  MODEL 

Let  us  consider  the  completely  general  continuous  review  model  for 
multi-item  B(Q, r )  Model  which  was: 

N_ 


minimize  Z=     V~         Wi         (5(ri) 


subject   to:  V~  .  ^i 


^         Ci    (  ri+     2    -  Ai)  4    K 


N 


£i    Ql 

Qi  >   0,    and  r^     unrestricted 


» 
where     ?  (  ri)  =  I  ['i2  +  (  ri-,fli)2]  §  (-H^-)  -  |  (  rl,*i)$  (^) , 

and        Wi  is  a  weighting  factor. 

A.   NECESSARY  CONDITIONS 

Applying  the  Lagrange  Multiplier  technique,  the  following  results  will 
be  obtained: 

n  W1   D  ri-  Qi    „       1       f^-'Ai 


L(Qi,    rif^,    9)    = 


J^Qi    P(ri>+^  Ci(ri+*2   -^i)-KJ    +6 K^  Qi  -^ 
i=l  lj-1  J         *i=l 

1L_  =-Wi_    A(ri)   +%i   .    M±  „  o    ;  (V.2) 

3Qi         Qi2    r  2  0^2 


(V.l) 


7i  ■  -f- [<  r^>  $  ( ¥r->  -«*=#>]    +4"  ■  o=      <v-3) 


^~  =  y"     Ci    (  ri  +  ^->fci)    -   Kn    =  0;    and 


(V.4) 
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N 


^-  =^%"   '  K~       '   °-  <V'5> 


-de     i—  qi        2 

i=l 


Then  equations  (V.l)  and  (V.2)  can  be  solved  for  each  multiplier  yielding 

2  ^i  Wi 


&  .   Qi^i+C^— >  "  <  rl-A)$  (^^)   -  2W^(  ri) 


or 


Qj  =    2CwiP(ri)  +8%i3 
^Ci 


^  .    ^(-s^fe-  (n-ioSHfe^: 


y        Ci(  ri  +  2  -JUL  i)  =  K   , 
1=1 


(V.6) 


(V.7) 


(V.4) 


and 

A 


IN 

z 


Considering  equation  (V.7),  the  right-hand  side  of  equation  is  always 
positive,  which  implies  /J^>  0,  since  4).  =  0  implies  Qi  =  00    which 
violates  equations  (V.5)  and  (V.4) . 

Observing  equation  (V.6),  the  right-hand  side  of  this  equation  is 
always  positive.   This  suggests  two  possible  cases  for  considerations. 

Case  I.   11^0,  0=  0  (Ignoring  the  reorder  workload  constraint). 
The  necessary  conditions  in  this  case  are: 

2     2  Wi  P(  ri) 
Qi  =    ^ci ,  <v-8> 
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St  §  c-^i)   -  (  ri  -Ai)  $  (_^ 


■Ai< 


%=  T   *1 *    *1     ,  (V.7) 

CiWiQi 


and       N 

1=1 


IN 

y  Ci  (  rA  +  ^  -  Ai)   =  Kt  .  (V.4) 


Case  II.    /B>  0,  Q}0  (Both  constraints  are  active). 
The  necessary  conditions  are  the  same  as  the  previous  page  with  Eq 
(V.6),  Eq  (V.7),  Eq  (V.4),  and  Eq  (V.5). 

B.   ITERATIVE  SCHEME  FOR  CASE  I 

Equations  (V.8)  and  (V.7)  can  be  thought  of  describing  two  curves  in 
the  Q-r  plane.   Equation  (V.8)  clearly  shows  that  order  quantity  approaches 
zero  with  increasing  reorder  point  to  infinity.   Furthermore,  from  equation 

<V-8)>    fr1"  ^  °    f°r  aU  i; 
i 


d( 
x  ^   2  4lciWip(TiT 


Taking  the  second  partial,  -* —  is  greater  than  zero  or  less  than  zero 

dri2 

depending  upon  the   sign  of  numerator,  which  is 

-2 
d2Qi      .        2^CiWi  B(ri)    Hr$[Mt)  -  [c/l4>(I^t)  -  (ri-Ai)i  (-^)J 


dtt2 


3 


pZ^Ci  Wi  £   (ri)  J     7 

From  equation   (V.7),   note   that   the  nonnegative  quantity 

[A*«^>-<rt.-*i>i«Sfc*i»] 

decreases  to  zero  as  r^  increases  to  infinity,  at  which  time  the  order 
quantity  will  be  zero.   However,  taking  the  first  and  second  partials 
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2 

from  equation  (V.7),    dri   /   0,   d  ri  ^  0  will  be  obtained; 

dQi  dQ2 


dri  =  .lciWi_     y   0       for  all  i, 
0Ti 


dQi      $(rJ 


and     d  ri  =  4lCiWi  <fi        S   0       for  all  i  . 

If  one  plots  the  two  curves,  Figure  2  will  be  obtained. 
By  setting  the  value  of  ri  to  zero  in  Eq  (V.7)  and  Eq  (V.8),  Qa  and 
Qg  will  be  obtained  respectively,  which  are 

*♦<-&>- <rl-#i>$<-^) 

Qa    =     — — 


4},CiWi 


and 


a       »[l<'iW|<- •&>  +  $*♦<- 4*4 


To  have  a  solution,  it  is  necessary  to  have  Qa^Qr*   Equating  Eq  (V.7) 
and  Eq  (V.8),  yields 

.2    2W1  P(ri)  .     [M  C-1^)  -  (n-M)  $  (-^r^)J  2 
1       %<*  (H,CiWi)2 

or 

[^cK^fa-c^i)^)]2  (v-9) 

2CiWiJ(5(ri) 

Notice  that  Eq  (V.9)  is  the  function  of  the  decision  variable  r  only.   A 
double  iteration  as  a  function  ofSfj ,  Q  and  r  is  required,  i.e,  iterate 
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Figure  2.   Q  as  Function  of  r  for  Case  I 
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Q,  r  for  a  given  value  of  the  multiplier  f(L   .   For  the  initial  value  of  /f]  , 

[<M>c-£>^<§<-£>]2 


S  =  min 


for  all  i,  pic 


k    -  ± 


20^   (r.=0) 


as  a  convenient  starting  point.   Then  solve 


equation  (V.9)  for  the  [r ,'s}  ,  compute  £q.'sj  using  Eq  (V.8),  and 

N         Q 
calculate  the  value  of  constraint  using  Eq  (V.4).   Let 


^Mr^-A^H. 


i=l 

A  binary  search  will  be  used.   If  h)k  ,  increase  ^  by  -£-   .   If  H  <  K  , 

1  4  1 

g 

decrease  ffi.  by  .   Compute  the  value  of  H.   If  the  increase  (or  decrease) 

4 

of/fl^has  not  caused  the  sense  of  inequality  sign,  increase  (or  decrease) 

/fV  by  the  same  amount  —r~    .   If  the  sign  of  the  inequality  has  changed,  then 

C 

reduce  the  increment  to  ,  and  increase  (or  decrease)  >/)  ,  solving  for  the 

8 

f  r  '  s  %    and  £  Q. 's  Jat  each  value  of  /ft   and  computing  the  value  of  H  until 

the  sense  of  inequality  changes  again.   Continue  until  H  =  K  or  until  H  is 
within  the  tolerable  region  of  K, .   This  method  will  converge  to  the  solution 
vector  rapidly. 

From  the  Kuhn- Tucker  theorem  [_3_]  ,  if  we  have  a  convex  objective 
function  and  a  convex  constraint  region,  the  necessary  conditions  are  also 
sufficient.   Since  the  first  constraint  under  consideration  is  linear  in 
r,  the  region  is  convex.   To  show  Z(r^)  is  convex,  it  is  required  that 
the  Hessian  is  greater  than  zero  or  equal  to  zero;  i.e., 


OZi 


*, 


-3^3  Qi 


2 


^0 


Let  0C(r±)   =    di$(   ri,"/tl) 
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•*L 


)  . 


Taking  first  and  second  partial  with  respect  to  each  decision  variable, 
the  following  will  be  obtained: 

3Zi  m  _   Ot  (rj)         _\_  m   _  P(rt) 
9ri      OT"  ,      Qi       Qi2 

&t      _  _i_  $  e~-)  -a2zi      .     2P^i) 

^r^rj    Qi   *    ai         ^Q^Qj       Qi3   , 

*£Zj  =  0((ri)  02Zi     B  d   (rt) 

^i^i  "    Qi2   •  ^    «Qi^  rt       Qi2    g 

From  Equation  (V.10), 


^r^rj^Q^Qj    ^-3^   -9^3  Qi     Qt        q^      q2   q2 

Hence  Equation  (V.ll)  is  greater  than  zero  or  less  than  zero  depending 
upon  the  sign  of  bracketed  quantity.   Under  these  considerations,  the 
convexity  of  Z^  is  unknown;  only  a  local  minimum  can  be  assured. 

C.   EXAMPLE  OF  THE  GENERAL  MODEL  FOR  CASE  I 

Once  again  consider  the  inventory  of  three  items  from  IV. C,  it  is 

assumed  that  the  distribution  of  lead  time  demand  is  normal  with  mean 

9  f- Vi 

/*£  and  variance  (ft      for  the  i   item. 


Item  1 

Item  2 

Item  3 

%i 

1000 

1500 

2000 

ci 

1 

10 

20 

At 

100 

200 

300 

*l 

100 

100 

200,  and  let 

Kx  =  $8,000  and 

K2  ■ 

=  15. 

34 


Reviewing  section  (V.A) ,  note  that  the  solution  vectors  must  satisfy 

2  _  2W.  ft(ri) 

1      %  Ct  ,  (V.8) 

Ma  ^4)(I^)"(ri-^i)J(I^rL)  (v.7) 


CiWiQi 


and 


N 
J~     ci(ri+-~  -jl^)  =  Kx  .  (V.4) 

i=l 

Using  the  search  scheme  described  in  section  (V.B),  this  example  was 
solved  with  the  results  shown  in  Table  1  utilizing  three  investment 
levels.   Comparison  between  the  General  Model  for  Case  I  and  Simplified 
Multi-item  Model  of  section  III  was  made  for  the  purpose  of  determining 
how  good  or  bad  the  simplified  model  performed  in  terms  of  time-weighted 
shortages  and  marginal  cost  imputed  to  the  first  multiplier.   These  results 
are  shown  in  Tables  2  and  3.   It  can  be  seen  that  the  Simplified  Model 
produced  time-weighted  shortages  which  were  from  25%  to  nearly  100%  larger 
than  the  General  Model  for  Case  I. 

D.   ITERATIVE  SCHEME  FOR  THE  GENERAL  MODEL 

Again  equation  (V.6)  and  V.7)  can  be  thought  of  as  describing  two 
curves  in  the  Q-r  plane.   Equation  (V.6)  shows  that  the  order  quantity 

approaches  a  fixed  value,  Q^  as  the  reorder  point  goes  to  infinity. 

dQi  , 
Furthermore,  from  equation  (V.6),  ^  0  for  all  i; 

dri 


dQ 
dr 


i/[^)-Mitf)] 
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^\^  ITEM 

variable\^ 

1 

2 

3 

REMARKS 

r 

340.6250 

342.1875 

506.2500 

Kx  =  8,000 

Q 

60.8684 

79.4553 

178.4906 

% 

0.0044 

0.0044 

0.0044 

Zi 

0.1343 

1.7529 

7.8755 

Z  =  9.7627 

r 

291.4063 

264.8438 

307.8125 

Kx  =  4,000 

Q 

69.4359 

101.1958 

247.2386 

u 

0.0154 

0.0154 

0.0154 

Zi 

0.5334 

7.7740 

37.9866 

Z  =  46.2940 

r 

400.0000 

425.0000 

687.5000 

Kx  =  12,000 

Q 

52.6381 

60.4982 

137.0932 

tl 

0.0007 

0.0007 

0.0007 

zi 

0.0194 

0.2229 

1.0100 

Z  =  1.2522 

Table  2.   The  General  Model  With  an   Inactive  Reorder 
Workload  Constraint  (Case  I) 
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\^  ITEM 
VARIABLE^. 

1 

2 

3 

REMARKS 

r 

234.3750 

264.0625 

453.1250 

Kx  =  8,000 
K2  =  15 

Q 

746.5022 

289.1190 

236.0648 

41 

0.0055 

0.0055 

0.0055 

z 

i 

0.2179 

2.7634 

10.5232 

Z  =  13.5045 

r 

175.7813 

181.2500 

296.8750 

Kx  =  4,000 
K2  =  15 

Q 

746.5022 

289.1190 

236.0648 

% 

0.0173 

0.0173 

0.0173 

zi 

0.8319 

11.4189 

42.8919 

Z  =  55.1428 

r 

300.0000 

343.7500 

612.5000 

Kx  =   12,000 
K2  =  15 

Q 

746.5022 

289.1190 

236.0648 

% 

0.0011 

0.0011 

0.0011 

h 

0.0386 

0.4632 

1.6465 

Z  =  2.1483 

Table   3.      Simplified  Multi-item  Model   Solutions 
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2 
(TQi 


The  second  partial,      y~   ,  is  greater  than  zero  or  less  than  zero 

dr* 

depending  on   the   sign  of   the  numerator; 


2  24t,CiWi    [P<ri.)  +   0fti]  $  (-§J**0  -  <*  (rt) 


2 


d  Qi  m  '»  I  ^v  i-/   w    J  ^  v  O'i 

dri  ["2/^CiWi  {  ^  (r  i)  +  G?-i  }  1  g 

Observing  equation  (V.7),  note  that  the  same  properties  as  section 
(V.B)will  be  obtained,  which  are 

£|  <  0,  and   &§   >  0   for  all  i. 

dQi  dQi2 

If  one  plots  the  two  curves  something  like  Figure  3  will  be  obtained. 

Observing  Figure  3,  the  numerical  search  in  Hadley  and  Whitin  [_l_j 
will  be  used  for  the  general  model.   To  initiate  the  numerical  procedure, 
it  is  necessary  to  determine  reasonable  multiplier  values.   These  two 

to   fl'l 

multipliers  are  used  to  compute  r  when  QL  =     u •■"   =  Q.   The  solution 
could  be  started  with  the  point  where  Q  =  Q  on  the  curve  defined  by  Eq 

JLi 

(V.6).   The  r  value  so  obtained  is  used  in  Eq  (V.7)  to  compute  a  new  Q 
value  Q2;  i.e.,  move  from  the  point  (Ql,  ri)  on  tne  curve  defined  by  Eq 
(V.6)  to  a  point  on  the  curve  defined  by  Eq  (V.7)  having  the  ordinate  r-^. 
The  Q2  value  is  used  in  Eq  (V.6)  to  compute  a  new  r;  i.e.,  we  move  from 
the  curve  defined  by  Eq  (V.7)  to  the  curve  defined  by  Eq  (V.6)  at  constant 
Q.  Hence  a  series  of  steps  is  obtained  as  shown  in  Figure  3.   It  is 
clear  that  the  numerical  method  must  converge  to  Q*  and  r*  for  every  i, 
given  QA^  QB  . 

Then  compute  the  value  of  Eq  (V.5)  and  Eq  (V.4). 
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r  * 


Eq  (V.7) 


Figure  3.   Q  as  Function  of  r  for  Case  II, 
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In  changing  the  values  of  two  multipliers  sequentially ,  ^  =  0.01, 
0.02,  ...  and  0  =  1,2,  ....a  set  of  solution  vectors,   {  Q,  rj  will 
be  obtained.   Then  compute  the  time-weighted  shortages  and  the  value  of 
Eq  (V.5)  and  Eq  (V.4).   'Good'  solution  vectors  can  be  found  with  small 
variance  of  Eq  (V.5)  and  Eq  (V.4)  and  approximate  the  minimum  time- 
weighted  shortages.   Once  obtaining  this,  iterate  with  more  careful 
selection  of  multipliers,  compute  the  objective  function  and  check  the 
variances  of  Eq  (V.5)  and  Eq  (V.4).   Continue  iterating  on  the  multipliers 
until  tolerable  differences  between  the  right  and  left  sides  of  Eq  (V.5) 
and  Eq  (V.4). 

More  sophisticated  approaches  are  available.  Among  these  is  a  search 

technique  proposed  by  Fiacco  and  McCormick  iji_\   ,  and  a  two  state  variables 

optimal  allocation  using  dynamic  programming. 

The  sequential  unconstrained  minimization  technique  [_4J  ,   SUMT,  is 

N 

E 


.  N 

based  on  the  minimization  of  a  new  function  P(X,f )  =  f (X)  +  pt""  l/gi(X) 


i=l 
over  a  strictly  monotonic  decreasing  sequence  of  P-  values   \'k  J  . 

Under  certain  restrictions  that  will  be  reviewed  subsequently,  there 

exists  a  sequence  of  feasible  points  4  X(fic)j   that  respectively  minimi 

{p(X  ,  p k)J  ,  and  it  follows  that  X  (  Pk)->  X  ,  a  solution  of  original 

function  as   v^  — >  0  (k-^OO  ).   The  following  is  a  concise  summary  of 

the  steps  describing  the  computational  algorithm: 

1.  Select  a  point  X   interior  to  the  feasible  region. 

2.  Select  f-p  the  initial  value  of  P  using 

Pj  =  -  Vf(x°)  T  v  p(x°)  /  |  vp(x°)|  2  , 


ze 
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or 


Pi- 


Vf(X°)  Hj  Vf(x°) 


-1 


"Vp(x°)V2  \Jp(x°) 


l2 

where    p(X)  =  ^J  .  I  /gi(X),  H-^  and  H2  are  the  Hessians  of  f(X)  and 
p(X)  respectively. 

3.   Determine  the  minimum  of  P(X,fj>)  for  the  current  value  of  Pk 
using  Gradient  Methods, 


x2  =  x    -     07p  (x1), 


or 

-1 


x2  =  x1  -  epaW)  ^x-^Xj]    y  pcx1) 

4.  If  k^|  ,  estimate  solution  using  extrapolation  formula. 

5.  Terminate  computations  if  final  convergence  criteria  are  satisfied, 
f(X)-  Gnt(£),  jU(t)J<  I   .   The  theoretical  optimum  value  vQ  is  bounded 
by  the  dual  and  primal  function  values  respectively, 

c[x(p),  JX(?)]£v0C      fCx(f)]  . 

If  the  value  v  is  not  within  the  bounds  above,  go  to  step  (6). 

6.  Select   V k+1  =  * k/C  >  wnere  c^  1« 

7.  If  k)  1,  estimate  minimum  for  reduced  p  -  value,  using  an  extra- 
polation formula. 

E.   EXAMPLE  OF  THE  GENERAL  MODEL  (CASE  II) 


Consider  an 

inventory  of  three  items  with  the 

following 

Item  1         Item  2 

Item  3 

U 

1000           1500 

2000 

Ci 

1             10 

20 

A 

100            200 

300 

<ft2 

100            100 

200 

Ki 

=  $8,000          K2  =  15 
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The  search  technique  proposed  by  Fiacco  and  McCormick  J_4j  was  used. 
The  following  initial  feasible  solution  was  used:   Q  =  600,  r,  =  200, 
Q2  =  270,  r2  -  260,  Q3  -  300,  r3  =  400.  21=^-^  .  lJMM   and 

*1  =  -  1.6991  were  obtained. 

This  example  was  solved  with  the  following  results, 


Qx  =  362.718,         r  =  322.1433, 
Q  =  266.403,         r2  =  272.0930, 

Q3  =  302.6870,        r3  =  425.7958,   P  =  -  0.00002231, 

and  the  value  of  the  objective  function  (Z)  was  Z  =  13.39781. 
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VI.   SUMMARY  AND  CONCLUSIONS 

It  is  apparent  that  the  general  models  proposed  do  not  represent  the 
ultimate  answer  in  multi-item  inventory  theory.   Since  the  conditions  of 
concavity  and  convex  sets  are  not  maintained  for  the  general  model,  the 
computational  algorithms  which  are  based  on  concavity  and  convex  sets  of 
objective  and  constraing  respectively  could  not  apply  directly.   Though 
the  P  function  in  section  (V.D)  appears  prohibitively  difficult  to  work 
with  computationally,  it  has  been  understood  that  it  is  minimized  efficiently 
and  accurately  for  the  great  preponderance  of  problems  solved  to  date,  by 
means  of  the  second  order  optimum  gradient  method. 

From  the  example  problem  we  see  the  general  problem  in  section  (V.E) 
provides  solution  which  is  better  than  the  solution  of  the  simplified 
continuous  review  model  in  section  (IV. C).   The  major  advantage  of  the 
simplified  continuous  review  model  is  its  computational  ease. 

While  an  efficient  algorithm  for  solution  of  the  general  problem  has 
not  been  presented,  the  advent  of  high  speed  computer  has  opened  this  field 
of  numerical  iterative  procedures  for  large  inventory  systems. 

Since  the  value  of  the  constraints  are  more  easily  determined  than  the 
order  cost  and  holding  cost,  the  model  proposed  seems  much  more  appropriate 
than  the  traditional  variable  cost  minimization  models. 

From  the  example  problem  given  in  section  (III) ,  it  is  apparent  that 
time-weighted  shortages  and  expected  number  of  units  short  per  unit  time 
seem  to  follow  each  other  closely. 

Finally,  all  the  solutions  to  example  problems  were  obtained  by 
utilizing  a  digital  computer. 
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